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Abstract. Using time-series of stereoscopic particle image velocimetry data, we 
study the response of a turbulent von Karman swirling flow to a continuous breaking 
of its forcing symmetry. Experiments are carried over a wide Reynolds number range, 
from laminar regime at Re — 10 2 to highly turbulent regime near Re — 10 6 . We 
show that the flow symmetry can be quantitatively characterized by two scalars, the 
global angular momentum I and the mixing layer altitude z s , which are shown to 
be statistically equivalent. Furthermore, we report that the flow response to small 
forcing dissymctry is linear, with a slope depending on the Reynolds number: this 
response coefficient increases non monotonically from small to large Reynolds number 
and presents a divergence at a critical Reynolds number Re c — 40 000 ± 5 000. This 
divergence coincides with a change in the statistical properties of the instantaneous 
flow symmetry I(t): its pdf changes from Gaussian to non-Gaussian with multiple 
maxima, revealing metastable non-symmetrical states. For symmetric forcing, a peak 
of fluctuations of I{t) is also observed at Re c : these fluctuations correspond to time- 
intcrmittencies between metastable states of the flow which, contrary to the very-long- 
time-averaged mean flow, spontaneously and dynamically break the system symmetry. 
We show that these observations can be interpreted in terms of divergence of the 
susceptibility to symmetry breaking, revealing the existence of a phase transition. 
An analogy with the ferromagnetic-paramagnetic transition in solid-state physics is 
presented and discussed. 

Keywords: Turbulence, Metastable states, Critical exponents and amplitudes 
(Experiments), Classical phase transitions (Experiments) 
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1. Introduction 

Symmetry breaking is an essential ingredient of classical phase transitions pQ. 
Remarkably, it also governs the transition to turbulence, that usually proceeds, as the 
Reynolds number Re increases, through a sequence of bifurcations breaking successively 
the various symmetries allowed by the Navier-Stokes equations coupled to the boundary 
conditions [2]. Finally, at large Reynolds number, when the fully developed turbulent 
regime is reached, all the broken symmetries are restored in a statistical sense [3]. 

Beside this well-established scenario, some experiments in closed flows have raised 
intriguing features proceeding during transition to turbulence [U El El CO El El QUI EH] ■ 
For instance, Tabeling et al. have evidenced a local peak in the flatness of velocity 
derivative in a von Karman flow of helium around Re = 2 x 10 5 [H] that was claimed to 
exhibit some characteristics of a second-order phase transition [9] and suggested to be 
associated with the breakdown of small-scale vortical structures. However, since then, 
no numerics or experiments ever confirmed this scenario, leaving this flatness peak as 
an unsolved problem. More recently, a transition has been evidenced in another von 
Karman flow around Re = 10 4 [TUl [TTj which consists in a global bifurcation of the 
mean flow from the basic symmetric flow topology towards two other flow topologies 
which spontaneously break the symmetry of the driving apparatus. This transition is 
associated with an hysteresis and a divergence of transition times, two features classically 
observed in phase transitions. Even more recently, in a von Karman flow very similar to 
the flow considered here, de la Torre and Burguete observed a time-intermittency of the 
large Reynolds number turbulent state, switching in between two metastable symmetry- 
breaking states [121 US]- Such transitions with dynamical symmetry breaking of mean 
patterns of the turbulent state are observed in 2D turbulence simulations (THUS], where 
they have been claimed to have strong geophysical and astrophysical relevance. They 
have also been observed and modeled in dynamo regimes of a liquid-sodium turbulent 

flow dam. 

Turbulent flows being intrinsically out-of-equilibrium systems, there is a priori 
no reason to describe them using tools borrowed from thermodynamics. However, in 
recent years, developments of out-of-equilibrium thermodynamics led, in particular, to 
a number of interesting applications to turbulence [IH1 EEl ED]. In this context, one 
may wonder whether the observed turbulent transitions can be interpreted in terms of 
phase transitions with a symmetry-breaking or susceptibility divergence signature. In 
the above mentioned works, while the Reynolds number appears as the natural control 
parameter, the equivalent of an order parameter and its associated susceptibility has 
not been clearly identified. 

In this paper, we consider a von Karman turbulent flow and its averaged large- 
scale patterns. We introduce a susceptibility to symmetry breaking for this flow and 
investigate its evolution as Re increases from 10 2 to 10 6 using Stereoscopic Particle 
Image Velocimetry (SPIV). We show that this susceptibility is as much as 50 times 
higher for the turbulent flow at very high Reynolds number than for the laminar flow. 
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Furthermore, in a intermediate Reynolds number range — 20 000 < Re < 200 000 — , we 
measure susceptibilities values as high as 300 times the laminar value and conclude to 
the existence, at a critical Reynolds number Re c = 40 000 ±5 000, of a divergence of this 
susceptibility to symmetry breaking. We also show that this phenomenon is associated 
with a peak in the amplitude of the fluctuations of the flow instantaneous symmetry. 
These observations complete the results of reference [21] with new complementary 
measurements at Re < Re c and with a re-processing of all the time-series. These new 
statistical analysis allow us to consider not only the mean values of the distributions 
but also their most-probable values. In particular, the divergence of the susceptibility 
and the peak of fluctuations which were previously found to occur at different values of 
Re [21] are now observed at the same Reynolds number if we consider the susceptibility 
based on most-probable values instead of mean values. This result sets the threshold 
for a "turbulent phase transition" at critical Reynolds number Re c = 40 000 ± 5 000. 
We relate the especially large fluctuations observed close to the transition to time- 
intermittencies associated with an ergodic exploration of a band of metastable states 
with spontaneous symmetry breaking. 

2. System, flow, symmetries and order parameters 

2.1. Experimental setup 

Our experimental setup consists of a Plexiglas cylinder of radius R = 100 mm filled up 
with either water or water-glycerol mixtures (cf. figure [lj a)). The fluid is mechanically 
stirred by a pair of coaxial impellers rotating in opposite sense. The impellers are flat 
disks of radius 0.925 R, fitted with 16 radial blades of height 0.2 R and curvature radius 
0.4625 R (cf. figure [lfb)). The disks inner surfaces are 1.8 R apart setting the axial 
distance between impellers from blades to blades to 1.4 R. The impellers rotate, with the 
convex face of the blades pushing the fluid forward, driven by two independent brushless 
1.8 kW motors. The rotation frequencies f\ and /2 can be varied independently from 1 
to 12 Hz, speed servo loop control ensuring a precision of 2%o for the mean frequency 
/ = (/i + /2)/2 and, for the relative frequency difference 9 — (fi — fzj/ifi + fi), an 
absolute precision of 1 x 10~ 3 and a stability along time better than 0.5 x 10~ 3 . 

Velocity measurements are performed with a SPIV system provided by DANTEC 
Dynamics. The cylinder is mounted inside a water or glycerol filled square Plexiglas 
container in order to reduce optical deformations. Two digital cameras are aiming at 
a meridian plane of the flow through two perpendicular faces of the square container 
providing 2D-maps of the three-components velocity field. Correlation calculations are 
performed on 32 x 32 pixels 2 windows with 50% overlap. As a result, each velocity is 
averaged on a 4.16 x 4.16 mm 2 window over the 1.5 mm laser sheet thickness. The 
spatial resolution is 2.08 mm. The data provide the radial u r , axial u z and azimuthal 
u v velocity components on a 95x66 points grid through time series of 400 to 27000 
fields regularly sampled, at frequencies from 1 to 15 Hz, depending on the turbulence 
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Figure 1. Schematic view of (a) the experimental setup and (b) the impellers blade 
profile. The arrow indicates the rotation sense, (c) 3D-view and symmetry : the 
system is symmetric with respect to any T^-rotation of angle 7r around any line in the 
equatorial plane which crosses the rotation axis. 

intensity and the related need for statistics. 

The control parameters of the studied von Karman flow are: 

(i) the Reynolds number Re = 2irfR 2 /i', where u is the fluid viscosity, which controls 
the intensity of turbulence and, 

(ii) the rotation number 9 = (fi — /2)/(/i + f'2), which controls the asymmetry of the 
forcing conditions. 

2.2. System symmetry vs. flow symmetry 

With exact counter-rotation of the impellers, i.e. when 9 = 0, the experimental system 
(cf. figure [IJa)) is symmetric with respect to any T^-rotation exchanging the two 
impellers: the problem conditions are invariant under 7r-rotation around any radial axis 
passing through the center of the cylinder (cf. figure pTc)). The symmetry group for 
such experimental system is 0(2) [22]. When the motors frequencies differ, i.e. when 
9 7^ 0, the experimental system is no longer 7?. T -symmetric and the symmetry switches 
to the SO(2) group of rotations. However, the parameter 9, when small but non-zero, 
can be considered as a measure of the distance to the exact 0(2) symmetry: the stricto 
sensu 50(2) system at small 9 can be considered as a slightly broken 0(2) system. Such 
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perturbative approach of the symmetry breaking has been successfully applied [231 El] 
for the 1:2 spatial resonance (or k — 2k interaction mechanism) with slightly broken 
reflection symmetry. 

Driven at a given 9 and a given Re, the system produces a flow. This flow — or 
at least its time-average at high Re — may be — or may not be — T^-symmetric. The 
problem is now to define a scalar quantity S, equivalent to an order parameter, which 
can quantify the distance of the flow to the 7^-symrnetry. The evolution of S with 9 
will measure the response of the flow to the symmetry of the forcing and x — 9S/ d9 will 
measure the corresponding susceptibility. Identification of relevant S can be made by 
studying the flow topology in the laminar and steady regime, at low Reynolds number, 
i.e. for Re < 200. 

2. 3. The von Kdrmdn flow topology: the laminar case 

When 9 = 0, the produced laminar flow is steady and composed of two toric recirculation 
cells separated by an azimuthal shear layer located at z = (cf. figure [2^ a)) reflecting 
the T^-symmetry of the system. Actually, each impeller pumps the fluid located near 
the rotation axis into its center before expelling it radially. This centrifugal pumping 
creates a toric recirculation cell in front of each impeller. In parallel, in each of the 
two toric cells, the fluid is rotating in the azimuthal direction following the nearest 
impeller. Then, when the motors frequencies differ, i.e. when 9 ^ 0, the shear layer 
moves towards the slowest impeller, breaking the symmetry of the flow, but keeping the 
two-cells topology (cf. figures [2tb-d)). As 9 tends to ±1, the flow continuously tends to 
a one-cell topology i. e. a topology with a single toric recirculation cell. 

2.4- Identification of two order parameters and associated susceptibilities 

We will use two different order parameters to quantify the distance of the flow to the 
7^-symrnetry. First, we consider a local quantity: the z-position of the shear layer of the 
flow, i.e. the position z s of the stagnation point on the rotation axi^fj The stagnation 
point position z s is only defined when the flow has the particular two-recirculation-cells 
axisymmetric structure and will not be suitable — because non-defined — to monitor 
the distance of the instantaneous flow to the T^-symmetry when it is turbulent. For 
this purpose, we prefer a second order parameter, a global quantity, the normalized and 
space-averaged angular momentum I(Re, 9, t) : 



where V is the volume of the flow. Practically, 1(f) is computed from SPIV- 
data restricted to a meridian plane only but, since azimuthal flow fluctuations are 

| Practically, for an axisymmetric flow, z s is measured as the axial position, ip(r — 0, z s ) = 0, of the 
zero iso-surface of the stream function ifj(r, z) which is defined through (u ri 0, u z ) = V x (r~ 1 ipe ip ) in 
cylindrical coordinates. 




(1) 
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Figure 2. Maps of velocity fields in the meridian plane (r, z) — the rotation axis is 
vertical at r = and the r -h- — r symmetry of the maps reveals the flow axisymmetry — 
for the laminar (Re = 1 20 ± 25) von Karman flow for different values of the symmetry 
control parameter 9: (a) 6 = 0, (b) 6 = -O.ff, (c) 6 = -0.19, (d) 6 = -1. The 
color maps, from blue to red ("jet" colormap), the azimuthal velocity u v , whereas the 
arrows map the (u r , u z ) field. The resolution of the fields has been reduced by a factor 
2 for better visibility. Positions r and z are given in units of vessel radius R. 



strong, time-average over several impeller rotation periods — statistically equivalent to 
spatial azimuthal averaging — estimate correctly the 3D-value of I(t) (see details in 
reference [25J). 

Examples of time-series of I(t) in turbulent regimes are provided below in figure [3] 
and reveal that strong spontaneous symmetry fluctuations can occur (cf. figure [3|a) for 
Re = 126 000, 9 = 0). We assume that ergodicity holds, meaning that the instantaneous 
turbulent flow is exploring along time its energy landscape according to its statistical 
probability. In this framework, the time average value I of I(t) is equivalent to a 
statistical mechanics ensemble average providing the average is performed over a long 
enough duration — comparable to the averaging time needed to measure z s — in order 
to correctly sample the slowest time-scales. The quantity / shows the major advantage 




Figure 3. Global angular momentum I(t) as a function of non-dimensional time 
f.t for experiments performed at Re = 126 700 ± 200 with (a) 6 = ± 10 and 
(b) 9 = -0.0084 and (c) at Re = 890 000 and 9 = 0. Green thin signals are SPIV 
data I(t) sampled at 15 Hz and blue thick signals correspond to low-pass filtered data 
If{t). Corresponding right figures shows the probability density functions (PDF) of 
the signals with same color code. The red lines indicate the most-probable values of 
I f (t). 



of being, contrary to z s , a robust observable since it is defined for each instantaneous 
velocity field whatever the structure of the flow and therefore even for turbulent flows. 
When observing the statistics of the fluctuations of /, we will sometimes encounter 
non-Gaussian and even multi-peaked distributions. Therefore, we will consider the time 
average J, the most probable value (s) J* and the fluctuations I(t) of the global angular 
momentum / to characterize the mean, most-probable and instantaneous symmetry of 
the flow. 

Then, using order parameter J, we define two susceptibilities of the flow to 
symmetry breaking xi an d Xi, based respectively on the mean and most-probable values 
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Figure 4. Experimental observation of proportionality between time-averaged global 
angular momentum / and mean shear layer altitude z s for (a) laminar flow Re = 
120 ± 25 (green squares, 7/z s = 0.235 ± 0.005) and (b) turbulent flows at Re = 67000 
(red circles, 7/z s = 0.227±0.005) and Re = 890 000 (blue triangles, 7/z s = 0.23±0.01). 
Each data series is obtained through plotting 1(8, Re) as a function of z s (9,Re) for 
varying 6 and constant Re. At high Re, this proportionality holds for —0.1 < 9 < 0.1, 
i.e. while the shear layer position z s is not saturated to the impellers positions ±0.7 R. 



of /, as: 



Xi 



m 
do 



and x~i 



de 



Another susceptibility may also be defined using z s as: 



Xz 



de 



(2) 



(3) 



Finally, since our observations (cf. figure |4j) have shown that / is strictly 
proportional to the mean altitude z s of the shear layer, with I / z s = 0.23±0.01 whatever 
Re, we may simply use the susceptibilities X = Xi = 0.23x 2 and x — Xi- m t ne sequel, 
we investigate the influence of turbulence on these order parameters — z s and I(t)— 
and the associated susceptibilities as Re, and therefore turbulence intensity, increases 
from 10 2 to 10 6 . 



2.5. About the statistics of the time-series and high-frequency noise 

The statistics of the long time series of I{t) obtained from SPIV data represent a crucial 
step in our study but they are limited by storage space and computing time. Our results 
are partly based on the mean values / and its equivalent z s computed on mean velocity 
fields and these lowest moments are well converged for all the data presented. Since 
statistical physics generally considers the most probable value of a distribution rather 
than its mean value and since some signals appear multivalued (cf. figure [31(a) and figure 



Susceptibility divergence, phase transition and multistability of a turbulent flow 



9 



gbelow), we are also interested in measuring absolute and relative most-probable values. 
However, the PDFs of the all raw signals I(t) are almost Gaussian and symmetric. Low- 
pass-filtering applied on I(t) with a typical 0.15Hz-cutoff removes some high-frequency 
noise, due in particular to the small scales of turbulence, and leads to PDFs with relevant 
non-Gaussian shapes. The best example is the high-resolution time-series corresponding 
to the highest fluctuation-level measured at Re = 39 000, presented in figure [8] The 
most-probable values are computed accordingly, but because of statistical limitations, 
these value are less converged than mean values. 

2.6. Formal analogy with ferromagnetic systems 

In the non-fluctuating laminar case, when 9 = 0, I is strictly equal to zero due to 
the symmetry of the flow. In contrast, as 9 drifts away from 0, the value of the 
angular momentum I becomes more and more remote from zero as the asymmetry 
of the flow grows. In the turbulent regimes, the same behavior is observed but with 
higher susceptibilities to symmetry breaking, i.e. higher sensitivity to 9. 

In such a framework, the Reynolds number — or a function of it [26] — is the 
equivalent of a statistical temperature and we are entitled to propose the following 
formal analogy between ferromagnetic systems and von Karman systems: 

• Order parameter: magnetization M angular momentum J; 

• Symmetry breaking parameter: external applied field h relative driving 
asymmetry 9; 

• Control parameter: temperature T <^ Reynolds number Re, or a function of it. 

3. Laminar and turbulent flow: influence of 9 and Re 

3.1. The laminar susceptibility 

The evolutions of the symmetry parameters z s and J as a function of 9 in the laminar 
flow at Re = 120 ± 25 are provided in figure [ij In this case I(t) = I(t) and z s are 
defined for each instantaneous field. We see that z s evolves almost linearly with 9 from 
z s = —0.7R (position of impeller 1) for 9 = — 1 to z s = 0.7 R (position of impeller 2) 
for 9 = 1, through z s = for 9 = (cf. figure |5]( a)). Variation of / with 9 is similar. 
In the central region, the linearity is excellent and we measure Xz — 1-03 ± 0.08 and 
Xi = 0.240 ± 0.005. 

3.2. Flow topology : the turbulent case 

Increasing the Reynolds number, one expects to reach fully developed turbulence around 
Re = 10 000 [UJ. In this turbulent regime and at 9 = 0, the T^-symmetry is broken 
for the instantaneous flow. However, as usually observed for classical turbulence, this 
symmetry is restored for the time-averaged flow (cf. figure [6^a)). Then, as in the case 
of the laminar flow, when 9 is varied, we observe the breaking of the 7^-symmetry of 
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Figure 5. (a) axial position z s of the stagnation point and (b) global angular 
momentum / = / as a function of 9 for the laminar flows at Re = 120 ± 25. 
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Figure 6. Maps of mean velocity fields of the turbulent von Karman flow at 
Re = 890 000 for different values of 9: (a) 9 = 0, (b) 9 = -0.0036, (c) 9 = -0.0147, 
(d) 9 = — 1, with same layout as figure [2j The r f-> — r symmetry of the maps reveals 
that the time-averaged mean fields are axisymmetric. 
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Figure 7. Turbulent regimes at Re = 67 000 ± 2 000 (red circles) and Re = 
890 000 ± 140 000 (blue triangles): axial position of the stagnation point z s (a-b) and 
global angular momentum / — mean I(t) (c-d) and standard error 07 (e-f) — as a 
function of 9. Wide span —0.3 < 8 < 0.3 on the left side (a,c,e) and xlO zoom 
—0.03 < 9 < 0.03 on the right side (b,d,f). Data for the laminar regime at Re = 120 
(green squares) from figure [5] are also plotted for comparison. The data have been 
symmetrized with respect to 8 = 0. Dashed lines in (d) indicates the tilted plateaus 
independent of Re and their extrapolations to ±7q- 



the turbulent mean flow, the structure of which matches quite well the laminar flow 
topologies, however with a higher sensitivity to 9 (cf. figures [6tb-d)). 



3.3. The turbulent susceptibility 

The main properties of z s and I with respect to 9 can be observed in figure [7] 
which combines the measurements obtained at an intermediate Reynolds number 
Re = 67000 ± 2 000 (red circles) as well as at the highest reachable Reynolds number 
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in our experiment Re = 890 000 ± 140 000 (blue triangles). The left-hand-side figures 
give a global view over a large span in 9: at this scale the two high- Reynolds- number 
series appear very close to each other and show much higher susceptibilities than the 
laminar flow series (green squares): turbulence enhances dramatically the sensitivity 
of the flow to symmetry breaking. The turbulent susceptibilities are at least 10 times 
larger than the laminar susceptibility since z s saturates — the shear layer disappears 
into one impeller — for \9\ larger than a finite and small value 9 C = 0.10 ± 0.02, close 
to the value predicted in [27] for flat disks and 10 times smaller than in the laminar 
case where the shear layer disappears asymptotically as 9 — > ±1 (cf. figure [2] and [5]( a)). 
These disappearances of the shear layer into impeller 1 or impeller 2 correspond to the 
transitions of the flow topology from two recirculation cells — at small \9\ — to a single 
recirculation cell — at large \9\. In figure [7(e), we observe two maxima for o~j, the level 
of fluctuations of J, around ±9 C , which are characteristic of these transitions. Such 
transitions have already been studied for other counterrotating von Karman flows in 
subcritical cases pm [25] . 

The effect of the Reynolds number on turbulent flow responses is revealed by 
the x 10- zooms in 9 on the right-hand-side of figure [7j It is clear that, in a narrow 
region \9\ < 0.01, we encounter two major differences between the Re = 67000 and 
Re = 890 000 experiments: 

• The susceptibilities are very different: xi — 42 ± 1 at Re = 67000 and xi — 9± 1 at 
Re = 890 000 (cf. figure [7(b,d)), respectively 175 and 37 times the laminar value. 

• The fluctuation level o~i(9) presents a sharp and narrow peak at 9 = for 
Re = 67000, but not for Re = 890 000 {cf. figure [7(f)). 

This peak is typical of all data in the intermediate Reynolds number range 20 000 < 
Re < 200 000 where the susceptibility is clearly higher than around the highest Reynolds 
number reached Re ~ 10 6 . Note that outside the active region of high susceptibility and 
high fluctuations \9\ < 0.01, there is no more .Re-dependency for 1(9): in the medium 
range 0.01 < \9\ < 0.05, not too close from the two-cells/one-cell transition at ±9 C , 
all curves collapse on tilted plateaus which — if extrapolated — cross the 9 = axis at 
I Q ~ ±0.05 whatever Re. 



3.4- Fluctuation level near 9 = 

For all Reynolds numbers in the intermediate range, we observe a sharp and narrow 
peak for the fluctuation level o~i(9) around 9 = 0. This is illustrated by two time-series 
at Re = 127000, in figure |3]( a) and (b), for 9 ~ and 9 = 0.0084 which correspond 
respectively to the top and the bottom of the fluctuation peak in figure [7(f). The 
difference between the two signals is striking and reveals that the fluctuation peak is 
due to the presence non-Gaussian intermittencies. These intermittencies have been 
observed in the range \9\ < 5 x 10~ 3 exclusively. The I(t) time-series at Re = 890 000 
and 9 = 0, in figure [3(c), corresponds however to the minimum fluctuation level of figure 
[7Tf): the corresponding PDF is narrow, symmetric and has a single central maximum. 
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Figure 8. (a) Global angular momentum I(t) for an experiment performed at 
Re = 39 000 for 9 = 0(t) = 0. The green thin line is SPIV data I(t) sampled at 
15 Hz and the blue thick line corresponds to 0.15 Hz low-pass filtered data I/(t). 
(b) Probability density functions (PDF) with the same color code. Blue dashed lines 
indicate the two most-probable values ±7* of If(t) and green dash-dot lines indicate 
the mean value I. (c) Time-series of 9(t) as produced by the motor regulation in unit 
of 10~ 3 : the correlation with I(t) is negligible (see text). 

3.5. Mean vs. most-probable values for highly intermittent regimes 

In the above paragraphs and figures we have presented the evolution of mean values of 
I(t) with 9 for different Re. Globally, most-probable values differ from mean values only 
when the susceptibility is high and are, in absolute value, greater than mean values. This 
difference is especially noticeable for Re around 40 000. The data set which corresponds 
to the maximum fluctuation level is presented in figure [8] where I ~ 0. It shows strong 
temporal intermittency between states of different angular momentum values revealed 
by the multi-peak PDF. 

The most-probable values of I(t) are plotted in figure [9j with respect to 9 very 
close to 9 = 0. We observe that for Re = 37000 ± 2 000, the slope x f° r most-probable 
values is approximately twice the slope x f° r mean values at 9 — 0. We also reported 
points obtained from three high-resolution time-series at Re = 39 000 (figure [sj) and 
43 000. For these data sets, we plot the first two most-probable values ±i* of I(t) 
corresponding to the two almost symmetric peaks of the PDFs with are exchanging their 
relative importance around 9 = 0. These values are significantly higher than those of 
the Re = 37000-data series. This can be due to either the proximity of the critical value 
Re c or to the limited statistics of the Re = 37000-data series. Anyhow, these results 
suggest that the corresponding susceptibility x is diverging at Re = 41 000 ± 2 000 as it 
will be discussed in the next section. 
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Figure 9. Mean / (red circles) and most-probable I* (blue stars) values of the global 
angular momentum I(t) for Re ~ 40 000. Close symbols come from a set of moderate- 
resolution data at Re = 37000±2000 (/ = 12Hz, 1 200 SPIV fields acquired at 1.7Hz), 
whereas open symbols come from high-resolution time-series at Re — 41 000 ± 2 000 
(/ = 8Hz, 9 000 or 27000 SPIV fields acquired at 15Hz as, e.g., in figure [8} acquired 
in a different glycerol-water mixture. The data have been symmetrized with respect 
to 6 = and, for the data point at ~ of the second series, we plot the two 
most-probable values ~ ±i* corresponding to the two peaks of the PDFs as in figure 



3. 6. Divergence of the susceptibility at Re ~ 40 000 

From the above results it is clear that the effect of the Reynolds number is to be studied 
for 9 = 0. The present section deals with the susceptibilities and the fluctuation level 
at 9 = and for Re covering the whole study range 10 2 — 10 6 . 

In figure [lO^a), we plot susceptibility measurements — x an d X, cf. equation [2] — as 



a function of Re. Figure 10 b) shows the inverses of these values and figure |l0[c) the 
corresponding standard deviation 07 of I(t). The susceptibility \ based on the most 
probable values I* of I(t) and the fluctuations of I(t) are both maximum for Re ~ 40 000. 
However, as reported earlier |2T] . the susceptibility x based on the mean values of I(t) 
reaches a maximum for higher Re ~ 90 000. We also note that x > Xi whatever Re, 
and that the highest measured susceptibilities are of the order of the highest measurable 
values considering the ^-precision of our setup. 

All the above results suggest a divergence of x(Re) which clearly appears with 
the fits in figures flOTa) and (b): x(Re) diverges at Re = Re c = 40 000 ± 5 000 with 



critical exponent —1 [illustrated in figure 10 with respect to l/log(i?e) (see discussion 



in section 6.2)] and resembles a classical magnetic susceptibility divergence at the 
ferromagnetic/paramagnetic phase transition. In the next section, we investigate further 
this transition by looking at the equivalent of the spontaneous magnetization of magnets, 
i.e. how, around Re c , the angular momentum I(t) spontaneously and dynamically 
transits between different finite values. 
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Figure 10. Reynolds number dependence of the symmetry-breaking susceptibility of 
the von Karman flow at 9 = (a) and its inverse (b). Black open symbols correspond 
to susceptibility \ — Xi — 0.23xz obtained from mean values I (squares) and z s 
(diamond); red full circles stands for \i obtained from most-probable values of 
/. The red full star corresponds to the region (Re = 41000 ± 2000) where infinite 
susceptibility is inferred from the strong multistability of the signal (cf. figures [8] and 
[9|. Dashed lines are fits of the most-probables values above and below the star as 
linear functions of l/log(i?e) as proposed in the discussion of section 6.2 (c) Standard 
deviation oj of the unfiltered global angular momentum I(t) at 9 = 0. 



4. Dynamic multistability near Re c and spontaneous "momentization" 

4-1- Observation of the dynamics at Re = 39 000 

Close to the susceptibility divergence, we observe a complex dynamics for the global 
angular momentum I(t), as illustrated in figure [8] at Re = 39 000 and 9 = 0. Indeed, 
one observes that I(t) does not just fluctuate randomly around zero — its mean value — 
but shows a tendency to lock for some times, preferentially around ±i* with of the 
order of 0.03-0.04, estimated from the main peaks of the PDF of the filtered signal If(t) 
(cf. figure §b)). 

In order to understand the nature of the fluctuations of I(t), we have first 
checked that I(t) is un correlated with the control parameter 9(t) (correlation coefficient 
C(I,9) = I(t) 9(t) joi oq < 0.02, cf. figures |8](a,c)) and thus reveals only spontaneous 
fluctuations. Actually, I(t) fluctuates with two separate time scales: (i) fast fluctuations 
related to "traditional" small-scale turbulence and (ii) intermittencies corresponding to 
residence time of typically few tens of seconds and certainly due to the proximity of the 
susceptibility divergence. If one performs a time-average of the velocity field over one 
of these intermittent periods, one obtains a time localized "mean" flow, analogous to 
what is obtained for true mean flows when 9 ^ as presented in figures [6]^b-c), and 
with its own level of spontaneous symmetry breaking. To go one step further with the 
correspondence between I(t) and the amount of symmetry breaking or the position of 
the shear layer, we have computed bin-averages of the velocity field for a time-series 
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Figure 11. Bin-averaging of the velocity field with respect to J(t) for Re = 39 000 
and 9 = (same data series as in figure [8]) . Three bin-averaged velocity fields are 
presented for bins: (a) -0.075 < 7/(t) < -0.050, (b) -0.0125 < I f (t) < 0.0125 
and (c) 0.0375 < If{t) < 0.0625. (d) Corresponding stagnation point axial position 
z s measured for 11 such bin-averaged fields for I between —0.075 and +0.075 (green 
stars), superimposed to regular time-averaged data for varying 8 at Re = 67 000 (red 
circles) and Re = 890 000 (blue triangles) as in figure |2]jb) . 



I(t) in 11 bins ranging from —0.075 to +0.075. In figures 11 a-c), we have plotted three 
examples of such bin-averaged fields. The pattern is ^-symmetric for the bin around 
zero (cf. figure [TTj^b)) and non-symmetric otherwise (cf. figures 11 'a) and (c)). From 
the resulting series of conditionally-averaged velocity fields we measure the stagnation 
point axial position z s and report it as a function of the mean bin value (cf. figure 
TTJd)). There is a very good matching between these bin- averaged data and the time- 
averaged data already presented in figure [4] for different 6: beyond the proportionality 
between z s and J, I(t) can thus definitively be used as a quantitative instantaneous 
symmetry measurement. For example, the two symmetrical most-probable values ±7* 
of I(t) correspond typically to \z s \ ~ 0.2 R. 

From the above observations, we can conclude that, very close to Re c , the 
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turbulent flow explores a band of metastable symmetry breaking patterns evidenced 
by —I* % If if) % +1*, where If if) is the 0.15Hz low-pass filtered value of I if), the three 
most- visited states being — 0, and +1*. Farther from Re c , signals and PDFs also 
reveal coexisting preferred values for I(t) when 9 is very small. 



4-2. The "momentization" 

The dynamical spontaneous symmetry breaking observed in our experiment can be 
seen as the coexistence of multistable/metastable steady states of different symmetry 
which are explored dynamically along time. By analogy with magnetization we propose 
to describe this as successive "momentization" patterns observed in time, thanks to 
ergodicity which arises from the very high level of fluctuations. 

With respect to magnetism, the main difference is that the spontaneous 
momentization occurs along time and mainly very close to the critical point. 
Furthermore, owing to the statistical limitations, we cannot conclude if momentization 
occurs only on one side of the transition — as it is for the ferro/para transition — or on 
both sides. 



4-3. Statistics of the dynamic multistability 

The metastable/multistable states of finite lifetime (r) observed in time-series such as 
in figure [8^ a) can be characterized by their angular momentum averaged over their 
lifetime and noted J r , to avoid confusion with the full-signal mean value I. The f 
are distributed along the angular momentum axis and we note X(Re, 9) the /-interval 
containing the distribution for experiments at {Re, 9). 

Multistability is systematically observed in the high susceptibility zone limited by 
Re G X Re ~ [20 000, 200 000] and 9 G X e ~ [-0.01, +0.01]. Whatever Re and 9 in X Re and 
Xq, we observe that X{Re, 9) is confined in [— J , +/o] (lo has been defined at the end of 
3.3 and is shown in figure [7]). When Re ~ Re c , the interval X(Re, 9) ~ [—1*, +1* 



section 



is symmetric and close to [— 1 , +I ] (cf. figure [8]): the whole band is almost populated. 
However, farther from Re c , the interval X(Re, 9) is smaller and evolves inside [— Iq, +Io] 
with 9. 



5. Complementary investigation about the azimuthal structure of the flow 

Before discussing our results in the next section, it is important to note that all the 
above results concern 2D-velocity-fields in a meridian plane: the azimuthal spatial 
structure is ignored by our 2D-SPIV process since the azimuthal patterns — moving 
and fluctuating — are averaged along time and not resolved in space at all. In this 
section, we investigate some preliminary experimental results about the 3D-structure of 
the flow, gathered by visual observation with bubble seeding and white light [TTJ [25] . 
With a high density of bubbles we mostly see the patterns on the outer cylinder which 
are the radial vortices due to the shear-layer destabilization. In practice, we experience 
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that, even in the flow is turbulent, the azimuthal number m of vortices can be a well 
defined integer over some ranges of Reynolds number and/or some periods of time. So, 
exactly as for small dynamical systems, changes in m can reveal an Eckhaus instability 
for the 3D-mean pattern [§} 

This dynamics has been observed for 22 000 < Re < 300 000. Despite the vortices 
are strongly fluctuating, we can report the following observations of the average number 
m of vortices around the perimeter of the shear layer at 6 = 0: 

• m = 3 for 22 000 < Re < 47000, 

• m = 4 for 86 000 < Re < 300 000. 

Between these two regions, the fluctuations are generally too strong to draw precise 
conclusions, but at least a clear observation of m = 4 has been made once at Re ~ 51 000. 
So, we conclude that some kind ofm = 3f>m = 4 Eckhaus transition occurs between 
Re = 47000 ~ Re c and Re = 86 000 ~ 2Re c , in the region of very high susceptibility. 

6. Discussion 

We have reported the experimental study, for a von Karman swirling flow, of (i) the 
response of the mean flow to a continuous breaking of the T^-symmetry of the system 
and (ii) the spontaneous symmetry fluctuations of the instantaneous flow, from the 
laminar regime at Re ~ 10 2 to the highly turbulent regime at Re ~ 10 6 . The divergence 
of the susceptibility x(Re) for Re = Re c = 40 000 ± 5 000 reveals the existence of a 
phase transition or a bifurcation. In the following, we discuss the implications of this 
phase transition for turbulence, the critical behaviors and the analogy with the para- 
ferromagnetic transition. 

6.1. A Reynolds- dependent turbulence 

This new transition makes some non-dimensional characteristic quantities of the flow still 
.Re-dependent at very high Re. The turbulence is thus definitively not fully developed in 
this closed turbulent flow, contrary to the most common observations in this Reynolds 
number range. In other closed flows, the literature reports some examples of transitions 
occurring at Reynolds number of the order of 10 5 , a decade in Re beyond where the 
turbulence already looks like fully developed. First, in various type of turbulent fluid 
experiments, several authors report transitions and symmetry breaking at high Reynolds 
number (cf., e.g., references jH [5j El [7]). Moreover, in a similar von Karman flow, de la 
Torre and Burguete [121 H3] observed a bistability between two broken 0(2)-symmetry 
states with z s /R = ±0.19. Although the slow dynamics is very different from what is 
reported here, both experiment are likely to belong to a common framework (see below). 

§ We use here the word "mean" not for the 2D- or axisymmetric-time-average, as everywhere in the 
present article: here, "3D-mean" involves time averaging in a rotating frame following the azimuthal 
pattern. Such image processing will be a challenge, implying movies of the whole system [28] or full 
3D velocimetry. 
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Finally, in a liquid helium von Karman flow, Tabeling et al. report a local peak in the 
flatness of velocity derivative around R\ = 700 — corresponding to Re = 2 x 10 5 [8]. 
An interpretation in terms of a second-order phase transition was proposed [9]. Velocity 
derivatives cannot be measured in our experiment because of the low spatial resolution 
of our SPIV system. However, if one could establish that the observations of Tabeling 
et al. are indeed connected with ours, then our work would prove that their transition 
is related to the global structure of the mean flow and thus does not result from the 
breakdown of small-scale vortical structures as proposed in reference [9]. 



6.2. Critical behaviors and analogy with magnetism 

6.2.1. The control parameter: a function of log Re - Since our study covers 4 orders 
of magnitude in Reynolds number, it is natural to acquire and plot results along a 
logarithmic scale. However, we could search for critical behaviors with either Re 
or log i?e as control parameter. Indeed, critical behaviors are asymptotic behaviors 
around a given threshold Re c and therefore both (Re — Re c ) and (logi?e — logi?e c ) 
are mathematically equivalent. However, it is very common in non-linear physics, 
that asymptotic behaviors or normal forms correctly describe the dynamics, not only 
asymptotically, but over finite — and even large — ranges around thresholds using one 
particular control parameter. This is exactly what happens here and because of the 



quality of the l/x(i?e) fit in figure 10 b), we conclude that a relevant control parameter 
for this critical phenomenon can be logi?e or a gentle function of it. 

Within the formal turbulence/magnetism analogy, we actually use Castaing's 
proposition [26J for a temperature of turbulent flows which is T ~ 1/ logi?e. 



6.2.2. Critical behavior of the susceptibility - The divergence appears clearly in figure 
[lO](b) where l/x(Re) goes to zero at Re = Re c = 40 000 ± 5 000 in a way that can be 
fitted with a —1 critical exponent by: 

1 1 

^ log Re log Re c 

with different prefactors below and above Re c . This expression can be interpreted as 
the combination of the Castaing statistical temperature of turbulence T ~ 1/ log Re 
with the classical mean-field-theory expression for susceptibility at a critical point: 

X oc |T - Tcf 1 , 

Expression ^ has already been successfully used in reference [21] to describe the 
"mean" susceptibility x computed with mean values — instead of most probable values 
for x — of I{t) and with a more limited set of data. As already underlined, this rather 
fair agreement between mean-field-theory and our turbulent measurements is at first 
surprising. However, recent theoretical developments suggest that in turbulent von 
Karman flow the interactions are long-range like [291 EQ] supporting therefore the idea 
that mean field theory is nearly applicable in our system. 



(4) 



Susceptibility divergence, phase transition and multistability of a turbulent flow 



20 



6.3. Momentization vs. magnetization 

Within the proposed analogy between our observations and the para-ferromagnetism 
transition, the equivalence between the magnetization and the momentization is the 
most questioning aspect. In particular, the shape of the J(#)-curves is more similar to 
the magnetization response curves of a single magnetic domains than to poly-domains. 
For the para-ferromagnetism transition, the magnetization increases with the 



distance below the critical temperature while here, as stated in section 4.3, the 
spontaneous momentization is finite, limited in amplitude by Iq, and decreases with the 
distance to Re c . The system may share some characteristics with reentrant noise-induced 
phase transition similar to that observed in the annealed Ising model [3T1 132] . Our 
system is definitively of finite size and it is thus coherent that both the level of fluctuation 
ai(Re c ) and the band of spontaneous fluctuations for I(t) — which corresponds to a finite 
excursion of the shear layer position along z — remain finite at critical point (figure 

Be)). 

In their von Karman experiment at 9 = and Re ~ 4 x 10 5 , de la Torre and 
Burguete ([121 [13]) also observe multistability of three and only three states with can be 
interpreted, within our framework, as f = —I*, and +/*. While the amplitude of the 
J T = mode is constant along time, their system switches between the two f = ±i* 
modes with a very slow dynamics (hours) and very fast transitions. This dynamics 
does not seem to depend on Re once Re >2x 10 4 , the symmetric I T = state being 
dominant below. Within our framework, we could propose that, in this experiment, the 
susceptibility is infinite whatever Re above 2 x 10 4 (cf. figure 5(a) of reference [12]). 



6.4- Specificity of turbulent fluctuations with respect to classical thermal noise 

In our turbulent system, the observed time-dynamics seems to be due to the very 
high level of intrinsic turbulent fluctuations, which are very different from classical 
thermal fluctuations that pilot the dynamics of magnetic momentums. These 
turbulent fluctuations are considered as the noise or the statistical temperature origin 
of our system. In classical close-to-equilibrium thermodynamical systems, thermal 
fluctuation are generally considered as regular additive noise exactly as for the classical 
treatment of instabilities in the presence of noise where an external noise — additive or 
multiplicative — is introduced as a perturbation of a regular system [33J. For turbulent 
flows in general and for our von Karman flow in particular, the situation is different: the 
flow itself is intrinsically highly fluctuating so that fluctuations are of the same order 
of magnitude than the mean values [25] and cannot be treated as perturbations. So, 
the specific problem our experiment addresses is the behavior of the instability of an 
average global quantity such as I(t) in an intrinsically fluctuating or noisy system. This 
is similar to the problem of the experimental dynamo instability [311 EH] which can only 
occur in turbulent flows because of the nature of the existing conducting fluids of the 
Universe. Therefore, non-classical behaviors are likely to be observed for bifurcations 
or phase-transitions in such highly-turbulent flows and our system represents a unique 
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tool to study these transitions. 
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